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Matrix polynomials, generalized Jacobians, 
and graphical zonotopes 

Anton Izosimov* 


Abstract 

A matrix polynomial is a polynomial in a complex variable A with coefficients in n x n 
complex matrices. The spectral curve of a matrix polynomial P{X) is the curve {(A,/i) G 
I det(P(A) — /r • Id) = 0}. The set of matrix polynomials with a given spectral curve C is 
known to be closely related to the Jacobian of C, provided that C is smooth. We extend this 
result to the case when C is an arbitrary nodal, possibly reducible, curve. In the latter case 
the set of matrix polynomials with spectral curve C turns out to be naturally stratified into 
smooth pieces, each one being an open subset in a certain generalized Jacobian. We give 
a description of this stratification in terms of purely combinatorial data and describe the 
adjacency of strata. We also make a conjecture on the relation between completely reducible 
matrix polynomials and the canonical compactified Jacobian defined by V. Alexeev. 
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1 Introduction 


The study of isospectral matrix polynomials lies at the crossroads of integrable systems and 
algebraic geometry. An n x n matrix polynomial of degree m is an expression of the form 
P(A) = Aq + AiX + • • • + AmA™ where Aq, ..., Am G Mat„(C) are complex nxn matrices, and A 
is a complex variable. It is well known that there is a close relation between matrix polynomials 
and line bundles over plane curves. With a given matrix polynomial P{X), one can associate 
a plane algebraic curve {(A,/i) G | det(P(A) — /r • Id) = 0}, called the spectral curve, and 
a line bundle over the spectral curve given by eigenvectors of P{X). Conversely, given a plane 
algebraic curve C and a line bundle E ^ C satisfying certain genericity assumptions, one can 
construct a matrix polynomial whose spectral spectral curve is C and whose eigenvector bundle 
is isomorphic to E. 

When the spectral curve is smooth, the relation between matrix polynomials and line bundles 
is described by the following classical result. Fix positive integers m and n, and let Me be the 
moduli space of n x n degree m matrix polynomials with spectral curve C, considered up to 
conjugation by constant matrices (for simplicity of notation, we omit the dependency of Me on 
m and n; note that the dimension n can always be found from the equation of the spectral curve 
C). Then, if the curve C is smooth and satishes certain conditions at infinity, the moduli space 
Me is a smooth variety biholomorphic to the complement of the theta divisor in the Jacobian 
of C. For hyperelliptic curves, this result is due to D. Mumford. In [18], Mumford used the space 
Me to give an explicit algebraic construction of the hyperelliptic Jacobian. For general smooth 
curves C, the space Me was studied by P. Van Moerbeke and D. Mumford in connection 
with periodic difference operators, and by M. Adler and P. Van Moerbeke |2|, A.G. Reyman 
and M.A. Semenov-Tian-Shansky |20|, and A. Deauville |9| in the context of finite-dimensional 
integrable systems. It is worth mentioning that the technique used in these papers is to a large 
extent based on earlier works of S.P. Novikov’s school on infinite-dimensional integrable systems, 
see, e.g., the review [12]. 

In the present paper, we are interested in generalizing the relation between isospectral matrix 
polynomials and Jacobian varieties to singular curves. It is known that for certain singular 
curves the space Me is a variety isomorphic to the complement of the theta divisor in the 
compactification of the generalized Jacobian of C, see e.g. m- However, in general, the space 
Me associated with a singular curve may even fail to be Hausdorff. For this reason, we replace 
the space Me with a slightly bigger space Pq defined as the set of n x n degree m matrix 
polynomials whose spectral curve is a given curve C and whose leading coefficient Am is equal 
to a hxed matrix B. The latter space is, by dehnition, an affine algebraic variety for any curve 
C and any matrix B (by an algebraic variety we always mean an algebraic set, i.e. we allow 
varieties to be reducible). Note that the space Me can be recovered from Vq by taking the 
quotient of the latter with respect to the conjugation action of the centralizer of B, provided 
that this quotient is well-defined. 

For a smooth curve C, the set was explicitly described by L. Gavrilov [13]. Namely, he 
showed that if C is smooth, then can be naturally identihed with an open dense subset 
in the generalized Jacobian of the curve obtained from C by identifying points at inhnity. We 
note that the latter Jacobian has a natural fiber bundle structure over the Jacobian of C, which 
agrees with the above description of the moduli space Me ss the quotient of . 

In this present paper, we consider the case of an arbitrary nodal and possibly reducible 
curve C. In this case, the isospectral set Vq turns out to have a rich combinatorial structure. 
This structure may be described in terms of orientations on the dual graph of C, or in terms of 
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a certain convex polytope, the so-called graphical zonotope m associated with the dual graph. 
Based on this combinatorial description, we make a conjecture on the relation of the set to 
the canonical compactihed Jacobian of C described by V. Alexeev [3]. 

The interest in varieties associated with singular curves also comes from integrable 
systems. There is a natural algebraically integrable system on matrix polynomials of the form 
Aq + Li\ + ■ ■ ■ + whose hbers are precisely the sets see [201 [g 13]. So,the 

results of the present paper can be regarded as a description of singular fibers for the integrable 
system on matrix polynomials. We also hope that these results will be useful on their own, in 
particular, for understanding the structure of compactihed Jacobians. 


Main results of the paper are Theorem 5.1 and Theorem 8.2 Theorem 5.1 describes the 
variety Vq associated with a nodal curve C as a disjoint union of smooth strata, each one being 


an open subset in a certain generalized Jacobian. Theorem 8.2 describes the adjacency of these 
strata. Two other results are Theorem 6.1 which describes the local structure of the variety Vq, 
and Theorem 


10.2 


which characterizes strata of Vq corresponding to irreducible and completely 


reducible matrix polynomials. Based on Theorem 10.2 we make a conjecture on the relation of 
the variety Vq to the canonical compactihed Jacobian of C (Conjecture 10.4). 


The main emphasis of the present paper is on constructions and examples. Detailed proofs 
will be published elsewhere. 


Acknowledgments. This work was partially supported by the Dynasty Foundation Schol¬ 
arship and an NSERC research grant. 


2 Statement of the problem 

Fix positive integers m, n, and let 

'Pm,n •= {P{^) = ^0 + + • • • + Afn-lX^ ^ + A^X^ \ Aq, . . . , Am G Mat„(C)} 

be the space of n x n degree m matrix polynomials. For each P G Pm,n, let 

Cp := {(A, /i) G I det(P(A) - /r • Id) = 0} 

be the zero set of its characteristic polynomial. The curve Cp is called the spectral curve 
associated with a matrix polynomial P. This curve can be regarded as the graph of the spectrum 
of P depending on the value of A. Further, hx a matrix B G Matn,(C) with distinct eigenvalues, 
and let 


Prn,n ~ {BiX) — Aq -|- AiA J- • • • J- Am-lX^ ^ + AmX^ G Pm,n \ Am — J?} 

be the space of matrix polynomials P G Vm,n with leading coefficient B. The isospectral variety 
corresponding to a plane curve C is defined as 

Pc ■■= {P e Pin I Cp = C}. 

The main problem of the present paper is the following: 

Problem 2.1. Describe the structure of the set for an arbitrary nodal, possibly reducible, 
curve C. 
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Figure 1: Newton polygon of a spectral curve 


The solution of Problem 2.1 is given by Theorems 5.1, 6.1 and 8.2 below. We begin with 
several preliminary remarks. First, note that if the matrices B and B' are similar, then the 
corresponding isospectral varieties and Vq are naturally isomorphic. For this reason, we 
may confine ourselves to the case of a diagonal matrix B (recall that eigenvalues of B are pairwise 
distinct). 

Further, note that the choice of m, n, and B imposes restrictions on curves which may arise 
as spectral curves. Indeed, let C be a curve given by Q{X, /r) = 0. Assume that is non-empty, 
and let 

P{X) = Ao + AiA + • • • + Am-iX^-^ + BX^ E vS- 
Introduce new variables z,w such that A = 1/z, and /r = vjjz^. Then 




z^O 


z z“ 


= lim ( det ( P - — • Id ) ) = 

2—>-0 


( 1 ) 


lim det(P + zA^-i + 

2-S>0 


+ z^^Ao -wld) = det(P - u; • Id). 


Geometrically, this condition means that the Newton polygon of the polynomial Q{X,fj,) lies 
inside the triangle with vertices (0,0), (0, n), {mn, 0), and that coefficients of the monomials lying 
on the hypothenuse of this triangle are equal to corresponding coefficients of the characteristic 
polynomial of the matrix B (see Figure]^. In what follows, we consider only those curves C 
which satisfy this condition. Note that for smooth and nodal curves, the above condition turns 
out to be not only necessary, but also a sufficient condition for a curve C to be the spectral 
curve for a suitable matrix polynomial P. Also note that for matrix polynomials of degree one, 
the above condition simply means that the curve C has degree n, and that its closure in CP^ 
intersects the line at inhnity at specific points prescribed by the eigenvalues of the matrix B. 

Finally, we note that there is a natural action of the group 

:= {U E PGL„(C) I UB = BU} 

on by conjugation. Provided that is non-empty, the quotient P^/G^ coincides with the 
moduli space 

Adc = {p E \Gp = G}/ GL„(C), 

where the group GL„(C) acts by conjugation. Indeed, there is a natural map tt; —)• AAc 

that sends every P E to its conjugacy class. Obviously, fibers of this map are precisely 
G^ orbits. Furthermore, since Pq is non-empty, there exists P E AAc such that its leading 
coefficient is conjugated to B. But this implies that this is so for any P £ AAc, because the 
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Figure 2: Nodal curves and their dual graphs 

eigenvalues of the leading coefficient of a matrix polynomial are uniquely determined by the 
spectral curve. Therefore, the map vr is surjective, and we have 

Me = V^/G^. 


Example 2.2. The following example shows that if the curve C is reducible, then the space 
Me is, generally speaking, not a variety. Let m = 1 , n = 2 , and let C be the union of two 
straight lines /r = ai + 5iA and /r = 02 + 62 A. To describe the space Me-, we first describe the 
variety Vq for B diagonal with diagonal entries 61,62 • A simple explicit computation shows 
that Vq is a disjoint union of three orbits Orbi U Orb 2 U Orbs where 


Orb2 


Orbi 


Q-i T 61 A 0 

0 02 + 62 A 


Ol + 61A 
0 


02 + 62A J j ’ 


Orbs 


oi T A 0 
z G C* 02 + 62 A 


Note that the closure of the orbit Orb 2 , as well as the closure of the orbit Orbs, contains the 
orbit Orbi, therefore the quotient Me = Vq/G^ is non-Hausdorff. 


3 Stratification of the isospectral variety 

Let G be an arbitrary nodal, possibly reducible, curve. In this section, with each matrix polyno¬ 
mial P G Pc’ associate a graph Tp and a divisor Dp on Tp, i.e. an integer-valued function 
on vertices of Tp. This construction leads to a stratification of Vq with strata indexed by pairs 
(a graph T, a divisor ZD on T) . 

Recall that dual graph Tc of a nodal curve C is defined as follows. The vertices of Tc are 
irreducible components of G. The edges of Tc are nodes. Two vertices are joined by an edge if 
there is a node joining the corresponding irreducible components. Nodes which lie in only one 
irreducible component correspond to loops in Tc (see e.g. Figure]^. 

Now we associate a subgraph Fp C Fc with each matrix polynomial in . Let P G Pe- 
Then, by definition of the spectral curve, for any point (A,^) G C we have dimKer (P(A) — p ■ 
Id) > 0. Moreover, if (A, p) is a smooth point of G, then dimKer (P(A) — p ■ Id) = 1, i.e. the p- 
eigenspace of the matrix P(A) is one-dimensional (see e.g. [6], Chapter 5.2). If, on the contrary, 
(A, ^) G C is a double point, then there are two possibilities: either dimKer {P{X) — p ■ Id) = 1, 
or dimKer (R(A) — p ■ Id) = 2 (cf. Section 2 of [T]). Now, define Fp C Fc as the graph obtained 
from the dual graph Fc by removing all edges of Fc corresponding to nodes {X, p) G G such 
that dimKer (R(A) — p ■ Id) = 2. Note that the subgraph Fp C Fc is generating, i.e. it has the 
same set of vertices as the graph Fc itself. 

Example 3.1. Let m = 1, n = 2, and let the curve C be the union of two lines Gi = {p = 
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ai + 6 iA} and 6*2 = {^ = 02 + 62 A}, as in Example 2.2 Let also i? be a diagonal matrix with 
diagonal entries bi, 62 - Then we have = Orbi U Orb 2 U Orbs (see Example 2.2). The only 


double point of the curve C is (Aq, /.io) = Ci H 6 * 2 , and the dual graph Tc is two vertices joined 
by an edge. For P G Orbi, we have dimKer (P(Ao) — /^o • Id) = 2, so the corresponding graph 
Tp is two disjoint vertices. For P G Orb 2 or P G Orbs, we have dimKer (P(Ao) — /io ■ Id) = 1, 
so Tp = Tc. 


As can be seen from Example 
generating subgraph Tp C Tc. 


3.1 


different components of may correspond to the same 


To distinguish between these components, with each matrix 
polynomial P G we associate an additional invariant, a divisor Dp on the graph Tp. This 
divisor is defined as follows. 

Let C'i,...,C'fc be the irreducible components of C. Consider the normalization X of C, 
i.e. the disjoint union of connected Riemann surfaces Aii,...,Afc, where Xi is the normal¬ 
ization of Ci- Note that X is equipped with two meromorphic functions A and /U coming 
from the embedding of the curve C into C^. For all but finitely many points u £ X, we 
have dimKer (P(A(u)) — n{u) ■ Id) = 1. This allows us to construct a densely defined map 
V^P: X —)• CP"'“^ which sends u to the one-dimensional subspace Ker (P(A(u)) — fi{u) ■ Id) C C"". 
Since X is smooth, the map ipp extends to a holomorphic map defined on the whole X. Now, 
consider a line bundle over X defined by taking the line 'ipp{u) for every point u £ X. Denote 
this line bundle by Ep, and let E*p be the dual bundle. The total degree of the bundle Ep can 
be found in the same way as in the smooth case (see e.g. [ 6 ], Chapter 5.2): 


^degEp Ixi = ^ (genus(Xi) - 1) + |E(rp)| +n 


( 2 ) 


2=1 


2 = 1 


where |E(rp)| stands for the number of edges of Tp, or, equivalently, the number of nodes 
(A, fi) £ C such that dimKer (P(A) — /r • Id) = I. For smooth curves, this formula reduces to the 
standard one 


deg Ep = genus(A) -|- n — 1. 


(3) 


Example 3.2. The following example shows that while the total degree of the bundle Ep can 
be computed in terms of the curve C and the graph Tp, the degree of the restriction of E*p to Xi 
can be different for different matrix polynomials P with the same spectral curve and the same 
graph. 

Let us again consider the two-lines spectral curve from Examples |2.2[ |3.I[ The corresponding 
Riemann surface is X = Xi U X2 where Xi is the closure of the line Cj = {/x = a* -|- biX} in CP^. 
Take a matrix polynomial P £ Orb 2 . Then for any point u = (A,/u) G Ci, we have 


P- H-ld 


0 2 

0 (02 -I- 62A) - (oi -I- 61A) 


so Ker (P — ■ Id) is spanned by the vector w = (1,0), which implies that the restriction of the 

bundle Ep to Xi is a trivial bundle, and degFp |xi= 0. Further, for u = (A,/x) G C 2 , we have 


P - /X • Id 


(ai-I-61A) - (a2-I-62A) z 
0 0 


so Ker (P — /x • Id) is spanned by the vector w = ( 2 ;, (02 -|- b 2 X) — (ai -|- 61 A)). The vector-function 
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0 


0 


fli “1“ 6iA 0 

0 02 + 62 A 


Ol + 61A * 

0 02 + 62 A 


Ol “t" b\X 0 
* 02 + 62 A 


> m -• 

0 1 


1 


0 


Figure 3: Graphs and divisors for matrix polynomials whose spectral curve is two straight lines 


w can be regarded as a meromorphic section of the bundle Ep over X2 with one pole of order 
one at infinity, so degEp \x2= “Ij and deg Ep |x 2 = 1- Analogously, for P G Orbs, we get 
degE*p |jVi= 1) and degE*p |x 2 = 0- Note that in both cases P G Orb 2 and P G Orbs, the total 
degree of the line bundle Ep is equal to one, as predicted by formula (§. 

Now, we use the numbers deg Ep \xi to define a divisor Dp on the graph Fp. Take any 
vertex Vi of Fp. By definition, this vertex corresponds to an irreducible component Ci of C or, 
equivalently, to a connected component Xi G X. Define 

Dp{vi) := deg E*p \x, - (genus(Ai) + n* - 1) (4) 


where n* is the number of poles of the function A on Xi or, equivalently, the degree of the 
equation of the irreducible component Ci in the variable Note that for smooth curves the 
right-hand side of Q is zero (cf. formula Q), so the divisor Dp measures the deviation of the 
degree of the line bundle Ep \xi from the degree that it would have in the smooth case. 

Example 3.3. Figure depicts the graph Fp and the divisor Dp for each of the components 
Orbi, Orb 2 , and Orbs of the variety Vq considered in Examples 


2.2 


3.1 


3.2 


Here and in what 


follows, dashed edges are those which do not belong to the subgraph Fp. 


Now, for any generating subgraph F of the dual graph of C, and for any divisor D on F, 
define the corresponding stratum 


5(F, D) := {P G I Fp = F, Dp = D], 


Example 3.4. When C is two lines (see Example 2.2 and Examples 
consists of three strata 


3.1 


B 


3.3), the variety 


( 0 -- *0 ) = Orbi, S ( 0 


1 ) = Orb 2 , S ( 1 


0 ) = Orbs. 


Theorem 5.1 below asserts that non-empty strata S{T,D) are smooth and connected, and, 
moreover, can be described as open subsets in semiabelian varieties. This result generalizes the 
description of Vq for smooth curves C found by Gavrilov [13j . The rest of the paper is devoted 
to the description of those divisors D for which S{T,D) is non-empty. We also describe the local 
structure of the variety in the neighborhood of each stratum and characterize those strata 
which are adjacent to each other. 
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4 Combinatorial interlude I: Indegree divisors and graphical zonotopes 

In what follows, a graph is an undirected multigraph possibly with loops. The notation C(r) 
stands for the vertex set of a graph T, and E(T) stands for the edge set. By Div(r) we denote 
the set of all divisors on T. Each divisor D G Div(r) may be regarded either as a formal integral 
linear combination of vertices of T, or as a Z-valued function on vertices. The degree \D\ of a 
divisor D is the sum of its values over all vertices. 

Let T be a graph, and let 0{r) be the set of all orientations on T. For each o € O, define a 
divisor 


indeg(o) := indego(u)u (5) 

ue v(r) 

where indeg|,(u) is the indegree of the vertex v in the directed graph (r,o), i.e. the number of 
edges pointing to v. 

Definition 4.1. Let T be a graph. We say that a divisor D G Div(r) is an indegree divisor if 
there exists an orientation o G 0{r) such that D = indeg(o). We denote the set of all indegree 
divisors by InDeg(r). 


Example 4.2. Figurej^depicts 8 possible orientations of a triangle and 7 corresponding indegree 
divisors. 


Remark 4.3. Divisors satisfying Definition 4.1 are also known in literature as score vectors m 
and ordered outdegree sequences [23] (note that there is no difference between indegree and 
outdegree divisors, since indeg(o) = outdeg(—o) where —o is the orientation obtained from o by 


inverting the direction of every edge). In algebraic geometry, divisors satisfying Definition 4.1 


are known as normalized semistable multidegrees. Such multidegrees were used by Beauville |8| 
to describe the theta divisor for a nodal reducible curve. The term normalized semistable 
multidegree is due to Alexeev [3]. Let us also mention a recent paper [5] where the authors 
consider so-called orientahle divisors. A divisor D G Div(r) is called orientable if there exists 
an orientation o G Cl(r) such that D = X^^(indegp(u) — l)v. Clearly, orientable divisors are in 
bijection with indegree divisors. 


The following result is well-known. 


Proposition 4.4. A divisor D on a graph T is an indegree divisor if and only if the degree of D 
is equal to the number of edges ofT, and for every subgraph Ti C T, the degree of the restriction 
of D to TI is greater or equal than the number of edges o/ Ti. 



1 


1 



1 


1 



1 


Figure 4: Eight orientations and seven indegree divisors on a triangle 
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For the proof of Proposition 4.4 see e.g. USKSIE]. See also [7] where it is shown that this 
result is equivalent to the well-known max-flow min-cut theorem in optimization theory. 

We will also need one more description of indegree divisors. To each vertex Vi of P, we 
associate a formal variable Xi = exp(ui), and consider the polynomial 


Br{xi, ...,Xn):= (Xi : 


[Vi,Vj\ 


where the product is taken over all edges of F. We have 


Br = P (exp(ui)exp(uj)) = ^ mult(r,T)) exp(i4), (6) 

[vi,Vj] DeDiv(r) 

where mult(r,D) is a non-negative integer for every divisor D. The following proposition is 
straightforward: 

Proposition 4.5. A divisor D = Y^diVi G Div(r) is an indegree divisor if and only if the 
polynomial B-p contains the monomial exp(D) = x^^ ■ ■ ■ x'jf, or, equivalently, if and only if 
mult(r,T)) > 0. Moreover, if D is an indegree divisor, then the number mult(r,i4) is equal to 
the number of orientations o G Cl(r) such that D = indeg(o). 

Example 4.6. Let P be a triangle. Then 


-Br = (xi -h X 2 ){xi + X3){X2 -k X 3 ) = 

x‘lx2 + Xixl + X^Xs + Xixl + X2X3 -I- X2xl + 2xiX2X3, 


( 7 ) 


which again shows that there are 7 indegree divisors on T (cf. Figure]^. 

Now, consider the space span]u(ui,..., Vk) = Div(r) (gjM formally spanned over M by vertices 
of F. Note that each edge of F may be viewed as a line segment in this space. The following 
definition is due to T. Zaslavsky, see [19]: 

Definition 4.7. Let F be a graph. Then the graphical zonotope Zp is the Minkowski sum of 
edges [vi,Vj\ G F considered as line segments in the space span]R(ui,... ,Vk)- 

Proposition 4.8. For any graph F, indegree divisors on F are exactly lattice points in the 
corresponding graphical zonotope Zp. 


Proof. It is well known that the Newton polytope of a product is the Minkowski sum of Newton 
polytopes of factors. Therefore, the graphical zonotope Zp is the Newton polytope of the 
polynomial Bp. In other words, Zp is the convex hull of all indegree divisors. In particular, any 
indegree divisor is a lattice point in Zp. Conversely, let B be a lattice point in Zp. Then D 
lies in the convex hull of indegree divisors, and, as follows from Proposition 4.4, it is itself an 
indegree divisor, q.e.d. □ 


Now we have four different descriptions of indegree divisors: in terms of orientations, in 
terms of linear inequalities, in terms of monomials in the polynomial Bp, and in terms of lattice 
points in the graphical zonotope Zp. 

Example 4.9. Let P = Kn be the complete graph on n vertices. Then 


Zp = Newton(xj -|- Xj) = Newton(xi — Xj) 


i<j 


i<j 


Newton(V(xi,..., Xn)) 
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( 0 , 1 , 2 ) ( 1 , 0 , 2 ) 


( 0 , 2 , 



Figure 5: Lattice points in the permutohedron P 3 


where Newton(i?) stands for the Newton polytope of a polynomial B, and V is the Vandermonde 
determinant 

V(xi,... ,Xn) = ]^ (xi - Xj) = det(x|“^) = ^ 

i<j a&Sn 


and Sn denotes all permutations of (0,1,..., n — 1). So, the graphical zonotope for the complete 
graph Kn is the convex hull of all permutations a € Sn- This polytope is known as the permu¬ 
tohedron Pn- Figure [^depicts the permutohedron P3 corresponding to the triangle K 3 . Lattice 
points in this permutohedron are seven indegree divisors (cf. Example 4.2 and Example |4.6[). 


Remark 4.10. Another way to define the graphical zonotope Zy is to consider the convex hull 
of indegree divisors of all acyclic orientations. Moreover, indegree divisors of the form indeg(o) 
where 0 is acyclic are precisely vertices of Zy (cf. |24) . Exercise 4.32). Note that if the graph 
F has loops, then it does not admit acyclic orientations. In this case, vertices of the graphical 
zonotope Zy are indegree divisors of those orientations that have no oriented cycles except 
compositions of loops. 


5 Description of strata 

In Section]^ we constructed a stratification = |J <S(F, D) where F is a generating subgraph 
of the dual graph Fc of C, and D is a divisor on F. Now, we describe those divisors D for which 
the stratum S{T,D) is non-empty. 

For a generating subgraph F of the dual graph Fc, denote by Cy the curve obtained from 
C by resolving those nodes which correspond to edges e G Fc \ F. Let also Cr / 00 be the curve 
obtained from Cy by identifying n points at infinity. 

Theorem 5.1. Assume that C is a nodal curve satisfying condition (§. Let also T be a gener¬ 
ating subgraph of the dual graph Fc of C, and let D G Div(F). Then 

i) the stratum S{T,D) is non-empty if and only if D is an indegree divisor on F; 

a) each non-empty stratum iS(r, D) is a smooth irreducible quasi-affine variety biholomorphie 
to an open dense subset in the generalized Jacobian of the eurve Cr / 00 . In particular, 

dim5(r,z)) = _ |i?(rc)| + |E(r)|; 

in) for each non-empty stratum iS(r, D), the quotient SiT, D)/G^ is biholomorphie to an open 
dense subset in the generalized Jacobian of the curve Cy- 
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( 0 , 1 , 2 ) - ( 1 , 0 , 2 ) 
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( 0 , 2 , 1 ) ( 2 , 0 , 1 ) 
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* 
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0 
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V3 / 
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[ * 

* 

yz J 


Figure 6: Strata corresponding to vertices of the permutohedron P3. 


Corollary 5.2. We have 

, mnin — 1) 

dimP^ =-^-1. 

2 

Remark 5.3. Details on generalized Jacobians of singular curves can be found in 12I1E2!. Here 
we recall that the generalized Jacobian Jac(C') of a nodal curve C is an extension of the Jacobian 
of the normalization of C with a commutative algebraic group (C*)^, where k is a non-negative 
integer. The same is trne for a nodal curve with identified points at infinity, provided that these 
points are all distinct. In particular, if C is a rational nodal curve, then Jac(C) ~ (C*)^, and 
Jac(C/oo) ~ (C*)™, where m> k. 


The proof of Theorem 5.1 is based on careful analysis of the correspondence between matrix 
polynomials and line bundles described in Section and is similar to Gavrilov’s proof in the 
smooth case [l 3 ]. Details of the proof will be published elsewhere. 


Example 5 . 4 . Let C be two lines, as in Example 2.2 and Examples | 3 . 1 |- 3.4 Then 


0*-*1 ) -<S{ 1*-*0 )-C*, 5 ( 0 -- •0)-a point. 


which means that for all strata we have an isomorphism S(T,D) ~ Jac(Cr/oo), i.e. in this 
example the open dense subset from Theorem | 5 . 1 | is the whole Jacobian. 


Example 5.5. Let m = 1, n = 3, and let C be the union of three straight lines Ci = = 

Qi + biX} in general position. Let also B = diag(6i, 62, bs). The dual graph Tc is a triangle, so 
|InDeg(rc)| = 7 (see Figures |^[^. Accordingly, the variety has seven strata of dimension 
three. Six of them correspond to vertices of the permutohedron P3; these six strata are depicted 
in Figure]^ The notation Ui stands for Oj + biX, and the numbers denoted by stars are assumed 
to be non-zero. Note that there is a one-to-one correspondence between these strata and Borel 
subalgebras of 013 that contain the Cartan subalgebra of diagonal matrices. 

Now, let us describe the seventh three-dimensional stratum 5 (rc, -I-U2 -I-U3) corresponding 

to the interior lattice point ( 1 , 1 , 1 ) in the permutohedron P3. Let 

(ci,C2,C3) := ( 1 , 1 , 1 ) X (61,62,63), k := (01,02,03) ■ (01,02,03) 


11 








1 


2 


0 


1 



Figure 7 : Four indegree divisors on two sides of a triangle 


where x is the cross product, and dot denotes the inner product. Then the stratum iS(rc, vi + 
V2 + Vs) consists of matrix polynomials of the form A + XB where an = Oj, off-diagonal entries 
of the matrix A satisfy 


012021 = C 3 Z, 013031 = C 2 Z, 023032 = ciz, 012023031 = w, 
and {z, w) / ( 0 , 0 ) is any point lying in the affine part of the nodal cubic 

w{kz — w) = ciC 2 Csz^. (8) 


Further, let us describe two-dimensional strata. There are three generating subgraphs of Fc 
with two edges, and each of them has four indegree divisors (see Figure]^, so has 3 X 4 = 12 
two-dimensional strata. These strata look similarly to strata in Figure]^ but with two stars (i.e., 
two non-zero elements) instead of three. Stars can be placed at any two off-diagonal positions 
which are not symmetric with respect to the main diagonal; there are exactly 12 ways to choose 
such two positions. Similarly, there are six one-dimensional strata that can be obtained by 
placing one star at any position, and one zero-dimensional stratum, the point diag(i/i, j^2) ^'3)- 


Remark 5 . 6 . Note that the stratum iS(rc, ui -|-112-1-113) in the previous example is isomorphic to 
(C*)^ times a thrice punctured sphere, which means that the open dense subset from Theorem 


5.1 


in this case is not the whole Jacobian. For any other stratum of Vq, we have an isomorphism 


S{T,D) ~ Jac(C'r/oo), as in Example 5.4 


Remark 5.7. Note that matrix polynomials lying in opposite vertices of the hexagon in Figurej^ 
are transposes of each other. More generally, if R is a symmetric matrix, then there is an 
involution a on given by cr{P) = P*. This involution maps each stratum 5(r, D) to another 
stratum 5(r, D') and thus induces an involution r: D ^ D' oxi the set InDeg(r). Explicitly, the 
involution r can be defined by the formula r(indeg(o)) = indeg(—0) where —0 is the orientation 
inverse to 0. In other words, we have t{D) = deg (F) — D, where deg (F) := X]j;eV"(r) 
the degree divisor of the graph F. Equivalently, the involution r can be described as the central 
symmetry of the graphical zonotope Z-p. 


6 Local structure and irreducible components 

In this section, we describe the local structure of the variety in the neighborhood of its 
smooth stratum S{r,D). Let a node be the germ at the origin of the complex-analytic variety 
{{z, w) e C‘^ \ zw = 0 }. 

Theorem 6.1. Let P G S{T,D). Then, in the neighborhood of the point P, the isospectral 
variety is locally isomorphic to the direct product 0/codimiS(F, R) = |R(Fc)| — |E'(r)| 
and a smooth disk of dimension dimiS(r,R) = dimP^ — |R(Fc)| -|- |R(r)|. 

^We use the notation codim 5(r, D) for the codimension of the stratum iS(r, D) in Pc, i-e- codim iS(r, D) ;= 
dim Vq — dimiSfr, D). 
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Example 6.2. Let C be two lines, as in Examples 2.2 and 3.1 - 3.4, and let 0 C Lc be a 


generating subgraph with no edges, i.e. the disjoint union of vertices of Lc. The only indegree 
divisor on the graph 0 is the zero divisor, and the corresponding stratum 5(0,0) is one point 
diag(ai + 6 iA, 02 + 62 A). Theorem 6.1 implies that near this point Vq is locally a node. Two 
disks forming this node are closures of the strata S ( 0 •-• 1 ) and 5 ( 1 •-• 0 ) • 


Example 6.3. Let C be three lines, as in Example 5.5 Then the stratum 5(0,0) is one point 


diag(z^i, 1 ^ 2 , t's) where = a* + 5jA. Theorem 6.1 implies that near this point the isospectral 
variety Vq is locally a product of three nodes, i.e. a union of eight three-dimensional disks. Six 
of these disks are closures of strata depicted in Figure and other two belong to the closure of 
the seventh stratum 5(rc', vi-\-V 2 + V 3 ) (note that this closure has a double point at the origin 
corresponding to a double point of the curve 


Corollary 6.4 (of Theorem 6.1). Irreducible components of the variety are in one-to-one 
correspondence with indegree divisors on the dual graph Tc. 


Proof. According to Theorem 6.1, each stratum 5(r,Zl) where T is a proper subgraph of Tc 
lies in the closure of a higher-dimensional stratum. Therefore, irreducible components of Vq 
are closures of strata of the form S{Tc, D), as desired. □ 

Example 6.5. Let m = 1, and let C be the union of n straight lines {/U = a* -|- ftjA} in general 
position. Let also B = diag( 6 i,..., 6 ^)- The dual graph Tc is the complete graph Kn on 
n vertices. Therefore, irreducible components of Vq are enumerated by lattice points in the 


permutohedron (see Example 4.9). It is easy to give an explicit description of components 
corresponding to vertices of the permutohedron: they generalize components depicted in Figure 
and correspond to n! Borel subalgebras of g[„ containing diagonal matrices. Explicit description 
of other components for n > 3 is unknown. 

Also note that the number of lattice points in the permutohedron is known to be equal to 
the number of forests on n labeled vertices. It is sequence A001858 in the online encyclopedia 
of integer sequences: 1, 2, 7, 38, 291,.... 


Remark 6.6. The idea of the proof of Theorem |6.1| is to show that if C is a nodal curve, 
then every P G Vq is a non-degenerate singular point for the integrable system on matrix 
polynomials (for the definition of non-degenerate singular points, see, e.g., [inj. Definition 7). 
Then one applies J. Vey’s normal form theorem |26j . 


7 Combinatorial interlude II: Multiplicities 


Definition 7.1. Let F be a graph, and let D be an indegree divisor on F. The multiplicity of 
D is the number mult(F,iA) entering formula ([^ for the polynomial B-p or, equivalently, the 
number of orientations 0 G OfT) such that D = indeg(o)[^ 

Example 7.2. Let F be a graph with 2 vertices joined by n edges. Then By = (xi + ^ 2 )’^, so 
multiplicities are given by binomial coefficients. 


Example 7.3. Let F = Kn be a complete graph on n vertices. Then 

V(xf,...,x2) 


Rr(xi, ...,Xn) = 


V(xi, . . . ,Xn)’ 


^Equivalence of these definitions follows from Proposition 


4.5 
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i.e. By is the Schur polynomial corresponding to the partition A = (n — 1, n — 2,..., 1, 0). 
Therefore, indegree divisors on T are weights of the irreducible representation of with high¬ 
est weight A, and the multiplicity of any indegree divisor is equal to the multiplicity of the 
corresponding weight. 


Remark 7.4. The multiplicity of an indegree divisor D can also be dehned as 1 plus the number 
of possibly disconnected oriented circuits in the directed graph (T, o) where o is any orientation 
of r such that indeg(o) = D, and a circuit is a cycle allowing repetitions of vertices but not 
edges. 

Note that this definition of multiplicities in terms of oriented circuits implies that for a 
loopless graph mult(r, D) = 1 if and only if D is a vertex of the zonotope Zy (cf. Remark 4.10). If 
the graph T has k loops, then for any D G InDeg(r), we have mult(r, D) > 2^, and mult(r, D) = 
2 ^ if and only if H is a vertex of the zonotope Zy- 


Now, let us define the multiplicity of an indegree divisor along an indegree divisor on a 
subgraph. Let Ti C T be a generating subgraph, and let ©(T \ Ti) be the set of all possible 
orientations of those edges of T which do not belong to Ti. Then, for each such partial orientation 
0 , we can define a divisor indeg(o) by the same formula ([^: the value of indeg(o) at the vertex 
V is the number of arrows pointing at v. Note that if Di G InDeg(ri), and o G 0{r \ Ti), then 
Di + indeg(o) G InDeg(r). 


Definition 7.5. The multiplicity of D G InDeg(r) along Di G InDeg(ri) is the number of 
orientations 0 G ©(T \ Ti) such that D — Di = indeg(o). We denote the multiplicity of D along 
Di as mult(r,L) | ri,iAi). 

Example 7.6. 



Corresponding partial orientations are 


/ ^^ 0 ' 


and 

T o' 


/ ;o' 


Remark 7.7. Note that for Ti = 0 and Di = 0, where 0 C T denotes a generating subgraph 
with no edges, Dehnition 7.5 reduces to Definition 7.1, i.e. mult(r,iA | 0,0) = mult(r,iA). 


Remark 7.8. Note that a necessary condition for the multiplicity of D along Di to be nonzero 
is D > Di- However, this condition is insufficient. For example. 



= 0 . 
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Remark 7.9. Another way to define multiplicities is by introducing a poset structure on the 
set of indegree divisors on generating subgraphs. Let F be a graph without loops. Consider 
the set IN(r) of pairs {T',D) where F' C F is a generating subgraph, and D G InDeg(F') is 
an indegree divisor on Fb Then the set IN(F) has a natural poset structure. Namely, say 
that (Fi,L)i) > (F 2 ,Z? 2 ) if Fi D F 2 , and there exists an orientation 0 G ©(Fi \ F 2 ) such that 
Di — D 2 = indeg(o). Let 'H(F) be the Hasse diagram of the poset IN(F). Then multiplicities 
can be defined by the formula 


mult(Fi,L)i I F 2 ,L> 2 ) 


number of directed paths from {T 2 ,D 2 ) to (Fi,L)i) in 'H(F) 
(|i?(Fi)|-|L;(F2)|)! 


It is worth mentioning that a closely related poset OVr was considered in m- The elements of 
OVr are pairs (F',D') where F' C F is a generating subgraph, and D' is a completely reducible 
indegree divisor on F' (see Section]^. The ordering is defined in the same way as above, i.e. 
OVr is a subposet of the poset IN(F). In [1], it is proved that the poset OVr is isomorphic to 
the facet poset of the Voronoi cell associated with the graph F. 


8 Adjacency of strata 


Now, we would like to understand which pairs of strata 5(F, D) are adjacent to each other. First 
note that an obvious necessary condition for the closure of the stratum iS(Fi,Z)i) to contain 
the stratum S{r 2 ,D 2 ) is dim5(Fi,L)i) > dim5(F2, £> 2 ), or, equivalently, |Ll(Fi)| > |i?(F 2 )|. 
Moreover, since for any double point (A, fi) £ C the set {L G V^ \ dimKer (T’(A) — /r • Id) = 2} 
is closed in V^, we should in fact have Fi D F 2 . However, the latter condition still does not 
imply that the closure of iS(Fi,L)i) contains iS(F 2 ,iA 2 ). 

To formulate a necessary and sufficient condition for a stratum to be contained in the closure 
of another stratum, let us first discuss in more detail the local structure of the variety V^ in 
the neighborhood of a stratum iS(F 2 ,iA 2 ). Theorem 6.1 implies that locally Vq can be written 
as 


(2C* U a point)^ x = |_| 2® 

1=0 


P 


{c*y X CL 


where p = codim5(F 2 , 112)5 and q = dimiS(F 2 , Zl 2 )- Thus, the neighborhood of the stratum 
iS(F 2 ,D 2 ) in the variety V^ is locally a union of 3^ smooth strata whose dimensions vary from 
q = dimiS(F 2 ,II 2 ) to p + q = dimV^. Each of these local strata is an open subset of some 
global stratum iS(Fi,Ili). 

Definition 8.1. The multiplicity of the stratum iS(Fi,Zli) along the stratum iS(F 2 ,iA 2 ) is the 
number of local strata S in the above stratification of the neighborhood of iS(F 2 , D 2 ) such that 
ScSiVuDi). 

Clearly, the multiplicity of iS(Fi,iAi) along iS(F 2 ,Il 2 ) does not depend on the choice of a 
point P G iS(F2,iA2). In particular, if this multiplicity is positive, then the stratum iS(F2,D2) 
lies in the closure of iS(Fi,Ili), and if the multiplicity is zero, then S{T 2 ,D 2 ) and the closure 
of 5(Fi,Di) do not intersect each other. Also note that if the multiplicity of iS(Fi,iAi) along 
S{T 2 ,D 2 ) is equal to one, then the closure of 5(Fi,L>i) is smooth at points of iS(F 2 ,iA 2 ). 
Multiplicity bigger than one means that the closure of iS(Fi,iAi) is singular along iS(F 2 ,D 2 ). 

The following theorem gives a necessary and sufficient condition for a stratum to be contained 
in the closure of another stratum, and also allows us to compute multiplicities. 
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Theorem 8.2. Let S{Ti,Di) and S{T 2 ,D 2 ) be two strata oJVq. Then 

i) S{r 2 ,D 2 ) lies in the closure of S{Ti, Di) if and only if T 2 C Ti and mult(ri,L)i 

r2, -D 2 ) > 0; 


a) if ^2 C Fi, then the multiplicity of the stratum 5(ri,F)i) along the stratum S{r 2 ,D 2 ) is 
equal to the multiplicity mult(ri,Di | r 2 ,-D 2 )- 


Example 8.3. Let C be three lines in general position, as in Examples |5.5| and |6.3[ Let 
also r = Fc, and let D = vi + V 2 + V'i ^ InDeg(F). We have mult(F,L)) = 2 (see 0 and 
Figure]^. Therefore, the stratum 5(F, D) has multiplicity 2 along the zero-dimensional stratum 
5(0, 0) = diag(ai + hiX). In other words, P = diag(aj + hiX) is a double point for the closure of 
SiT^D) (cf. Example |6.3| ). Moreover, for any (F',L)') / (0,0), we have mult(F,L) | T\D') < 1, 
which implies that the closure of S{T,D) is smooth at all points except diag(oi + hiX). 


Example 8.4. More generally, let C be a union of n straight lines (see Examples |6.5| and 7.3). 
Then maximal-dimension strata 5(F,L)) of the isospectral variety are indexed by weights 
of the representation of with highest weight A = (n — 1, n — 2,..., 1, 0), and the multiplicity 
of a stratum S{T,D) at the point diag(ai -|- biX) is equal to the multiplicity of the weight D. 
Also note that, in contrast to the case n = 3, some of the strata 5(F, D) may have multiplicities 
greater than 1 along strata of positive dimension, see, e.g.. Example |7.6[ 


Remark 8.5. Note that for any (Fi,iAi) and (F 2 ,D 2 ), one has mult(Fi,iAi | T 2 ,D 2 ) < 
mult(Fi,Di). Therefore, the closure of a stratum 5(Fi,iAi) is smooth if and only if 
mult(Fi,iAi) = 1, or, equivalently, if Fi has no loops, and Di is a vertex of the graphical 


zonotope (see Remark 7.4). 


Remark 8.6. Also note that conditions F 2 C Fi and D 2 < Di are necessary but not sufficient for 
the stratum 5(F2, D 2 ) to belong to the closure of the stratum 5(Fi, Di), see, e.g.. Example 7. 


9 Combinatorial interlude III: Reducible indegree divisors 

Let F be a graph. Recall that an orientation 0 G 0{T) is called strongly connected if any two 
vertices of the directed graph (F, 0 ) can be joined by a directed path; an orientation 0 G 0{T) is 
called totally cyclic if every edge of (F, 0 ) belongs to a directed cycle. It is well-known that if F 
is connected, then an orientation 0 G C(F) is totally cyclic if and only if it is strongly connected. 
On a disconnected graph F, an orientation 0 G C(F) is totally cyclic if and only if its restriction 
to every connected component of F is strongly connected. 

Proposition 9.1. Let D G InDeg(F) be an indegree divisor. Then the following conditions are 
equivalent: 

i) there exists no proper non-empty subgraph F' C F such that the restriction of D to F' is 
an indegree divisor; 

a) for every proper non-empty subgraph F' C F, the degree of the restriction of D to F' is 
strictly bigger than the number of edges ofV; 

in) there exists a strongly connected orientation 0 G 0 (F) such that D = indeg(o); 

iv) the graph F is connected, and D is an interior of the graphical zonotope Zr- 

®Note that if Zr is just one point, then we regard this point as an interior one. 


16 






Definition 9.2. Let D G InDeg(r) be an indegree divisor. If D satisfies conditions a)-d) listed 
in Proposition 9.1, then D is called irreducible. Otherwise, D is called reducible. 

Proposition 9.3. Let D G InDeg(r) be an indegree divisor. Then the following conditions are 
equivalent: 

i) the restriction of D to every connected component ofT is an irreducible dimsojQ' 

a) there exists a totally cyclic orientation o G 0{T) such that D = indeg(o); 

Hi) D is an interior point of the graphical zonotope Z^. 

Definition 9.4. Let D G InDeg(r) be an indegree divisor. If D satisfies conditions a)-c) listed 
in Proposition 9.3, then D is called completely reducible. 


10 Completely reducible matrix polynomials and compactified Jacobians 

Definition 10.1. A matrix polynomial P(A) is called reducible if there exists a proper non¬ 
empty subspace W C C"" which is an invariant under the action of the matrix P(A) for every 
A G C. A matrix polynomial P(A) is called completely reducible if every such invariant subspace 
W C C” admits an invariant complement. 

Note that if P(A) is reducible, then so is its spectral curve. However, the converse is in general 
not true. For example, matrix polynomials which belong to the stratum iS(rc,ui + V 2 + V 3 ) in 
Example |5.5| are irreducible, though the corresponding spectral curve is a union of three lines. 
It turns out that a necessary and sufficient condition for reducibility of P can be given in terms 
of the associated graph Pp and the divisor Dp: 

Theorem 10.2. Let P G <S(r, D) be a matrix polynomial. Then P is reducible if and only if D 
is a reducible divisor on P. Similarly, P is completely reducible if and only if D is completely 
reducible. 

Let us introduce the set 


{'Pc)cr ■= {-P(A) G Vc I P’(A) completely reducible}. 




By Theorem 10.2 we have 


{PE)cr = [_\^^^S(T,D) 


(9) 


where P is a generating subgraph of the dual graph Pc, and D is a completely reducible indegree 
divisor on P. 


Example 10.3. Let C be three lines, as in Example |5.5[ Then 



u 5 


Vo' 


cf. the explicit description in Example 5.5 


^Note that a restriction of an indegree divisor to a connected component is automatically an indegree divisor. 
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Now, we compare stratification © with the stratification of the canonical compactified Jaco¬ 
bian introduced by Alexeev. In [3] Alexeev showed that if C is a nodal, possibly reducible curve, 
then in degree 5 — 1 , where g is the arithmetic genus of C, there exists a canonical way of com- 
pactifying the generalized Jacobian. This canonical compactified Jacobian has a stratification 
that can be written in our terms as 


Jacg_i(C') = Jac(C'r) 


( 10 ) 


where T is a generating subgraph of the dual graph Tc, D is a completely reducible indegree 
divisor on T, and, as above, Cp is the curve obtained from C by resolving all nodes which 
correspond to edges not in T. 

Thus, from the combinatorial point of view, stratifications Q and (10) coincide. However, 
the strata themselves are different: each stratum of ([^ is biholomorphic to an open subset in 
the Jacobian Jac(Cr / 00 ), while strata of ( 10 ) are Jacobians Jac(C'r). To pass from Jac(C'r / 00 ) 
to Jac(C'r), we should take the quotient of S{T,D) with respect to the conjugation action of 
the centralizer of B (see Theorem 5.1). This observation leads to the following conjecture: 


Conjecture 10.4. Assume that C is a nodal curve satisfying condition ([^. Then the moduli 
space 


{'Pc)ctIG^ = {P € Vm,n I P is completely reducible, Cp = C} / GL„(C) 


of completely reducible n x n degree m matrix polynomials with spectral curve C, considered 
up to conjugation by constant matrices, is a variety isomorphic to an open dense subset in the 
canonical compactified Jacobian of C. 
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